Role of cross-helicity in magnetohydrodynamic turbulence 
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Strong incompressible three-dimensional magnetohydrodynamic turbulence is investigated by 
means of high resolution direct numerical simulations. The simulations show that the configuration 
space is characterized by regions of positive and negative cross-helicity, corresponding to highly 
aligned or anti-aligned velocity and magnetic field fluctuations, even when the average cross-helicity 
is zero. To elucidate the role of cross-helicity, the spectra and structure of turbulence are obtained 
in 'imbalanced' regions where cross-helicity is non-zero. When averaged over regions of positive 
and negative cross-helicity, the result is consistent with the simulations of balanced turbulence. An 
analytical explanation for the obtained results is proposed. 

PACS numbers: 52.35.Ra 



Introduction. — Magnetohydrodynamic (MHD) turbu- 
lence has been a starting point for modeling large-scale 
turbulent motion of plasmas in a wide variety of settings, 
ranging from laboratory experiments to astrophysical 
systems, [e.g., [lj]. When written in terms of the Elsasser 
variables the incompressible MHD equations read 
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where the Elsasser variables are defined as z^ — v ± b, 
v is the fluctuating plasma velocity, b is the fluctuating 
magnetic field normalized by y/Anpo, va = Bo/v^7rpo 
is the contribution of the uniform magnetic field Bo, 
P = (p/po + b 2 /2) includes the plasma pressure p and 
the magnetic pressure, and po is the background plasma 
density that we assume constant. It follows from these 
equations that for z =F (x, t) — 0, an arbitrary function 
z ± (x, t) = F ± (x± v^t) is an exact nonlinear solution 
that represents a non-dispersive wave propagating along 
the direction ^Fv^. Nonlinear interactions are thus the 
result of collisions between counter-propagating Alfven 
waves. 

The first phcnomcnological picture of MHD turbulence 
was proposed independently by Iroshnikov [2|, |3( and 
Kraichnan 4] (IK), who predicted an inertial range scal- 
ing of the isotropic energy spectrum E(k) ~ /c~ 3 / 2 . In 
this picture, the spectral energy transfer at a given scale 
A ~ \/k results from the cumulative effect of multiple 
weak interactions between counter-propagating Alfven 
wave packets moving along the magnetic field of the 
large-scale energy containing eddies. One shortcoming of 
this phenomenology is that it is based on the assumption 
of an isotropic spectral transfer, in clear contradiction 
with recent results that reveal the anisotropic charac- 
ter of MHD turbulence [e.g., [0, H, @- Indeed, Galtier 
et.al. 7] applied the formalism of weak turbulence to 
equations ([I]), demonstrating that the spectral transfer 
is much more efficient in the field-perpendicular plane, 
and derived a steeper energy spectrum E(k±) oc fcj" , 



where k± is the field perpendicular wave- number. This 
scaling was originally predicted in [8|, |9|, [lOj] based on 
more phenomenological grounds. 

As the cascade proceeds to smaller scales, the ed- 
dies become progressively more elongated in the field- 
parallel direction, and the nonlinear interaction becomes 
stronger. Eventually, the so called 'critical balance' 
condition of Goldreich and Sridhar (GS) [llj is estab- 
lished. This condition states that the turbulence is con- 
sidered strong when there is a formal balance between 
the crossing time of two interacting Alfven wave pack- 
ets and the characteristic nonlinear interaction time, i.e., 
kuBo ~ k±b\, where feii is the typical field-parallel wave- 
number of the fluctuations spectrum, and b\ is the rms 
magnitude of the fluctuations at the scale A ~ l/k±. The 
resulting scaling in the GS picture is E(k±) ~ kj_ . 

The explosive growth of massively parallel computers 
in recent years has made direct numerical simulations of 
MHD equations a valuable tool for studying fundamen- 
tal properties of MHD turbulence. For instance, sim- 
ulations indicate that the scaling of the energy spec- 
trum of strong MHD turbulence is anisotropic, as in 
GS picture, but with the scaling of the IK phenomenol- 
ogy 13, LL3|, LLJ, LUj, lift . In order to resolve this contro- 
versy, it has been proposed 17|, LL8| that magnetic and 
velocity fluctuations tend to align their polarizations in 
turbulent eddies, leading to scale-dependent depletion of 
nonlinear interaction of Alfven modes. This leads to the 



anisotropic energy spectrum E(k±) ~ k ± , in agree- 
ment with numerical simulations. 

The phenomenon of scale-dependent dynamic align- 
ment is closely related to the conservation of cross- 
helicity, an ideal invariant cascading toward small scales 
in a turbulent state. Cross-helicity has only recently be- 
come an object of systematic study, as it become clear 
that it plays a fundamental role in driven MHD turbu- 
lence 0, m Lnl, 0, d, HI HE LH LH- 

Denote E^ 1 = (|z ± | 2 )/4 the energy associated with the 
± waves, then the total energy and cross-helicity of the 



system are E = E + + E~ and H c = E + — E~ , respec- 
tively. Both energy and cross-helicity are invariants of 
the ideal MHD equations. Cross-helicity provides a mea- 
sure of the imbalance between non-linearly interacting 
waves. When it does not vanish the turbulence is called 
imbalanced, otherwise it is balanced. 

A significant interest to imbalanced MHD turbulence 
has also been motivated by astrophysical solar wind data, 
which indicate that solar wind turbulence is dominated 
by Alfven waves propagating outward from the sun [e.g., 
24| . A number of analytic derivations of the spectra 
of imbalanced MHD turbulence have been recently pro- 
posed 



Run 


Resolution 


Re 


a 


h/L± 


"1 


Pc 


A 


1024 2 x 256 


4000 


0. 


6 


1 





B 


512 2 x 256 


1500 


0.3 


5 


2 


0.6 


CI 


512 3 


1500 


0.25 


10 


2 


0.6 


C2 


1024 2 x 256 


4000 


0.25 


10 


2 


0.6 



TABLE I: Summary of simulations. A, B, C correspond to 
strong turbulence with different amount of net cross-helicity. 



dictory results. 
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[22j, which however lead to contra- 
Most of numerical simulations of MHD 
turbulence have so far concentrated on balanced cases, 
and practically no systematic study of strong imbalanced 
MHD turbulence in high-resolution direct numerical sim- 
ulations has been available. 

In an attempt to address the issue and to resolve the 
contradictions, we performed high resolution numerical 
simulations to investigate the inertial range of strong 
MHD turbulence with and without cross-helicity. Based 
on our results, we propose that in the imbalanced case the 
Elsasser energy spectra have different amplitudes, never- 
theless, their scaling is the same, E + (k±) ex E~(k±) oc 

3/2 

k ± . This scaling coincides with the scaling in the bal- 
anced case, which is consistent with the view that bal- 
anced MHD turbulence is as a superposition of locally 
imbalanced regions. We demonstrate that this picture 
is essentially consistent with the phenomenon of scale- 
dependent dynamic alignment, and provide an analytic 
explanation for the obtained spectra. 

Model equations. — The universal properties of weak 
and strong turbulence in MHD are accurately described 
by neglecting the parallel component of the fluctuating 
fields, associated with the pseudo- Alfven mode 16ll25|. 
By setting z,, = in equation (Q} we obtain the closed 
system of equations 



d_ 

di 



T v A ■ V|| z ± + (z* • Vjl) z ± = -VxP 



use this system to study strong anisotropic MHD turbu- 
lence. 

Numerical method. — We employ a fully dealiased 
Fourier pseudo-spectral method to solve equations ^ 
with a strong guide field (vA/v rms ~ 5) in a rectangu- 
lar periodic box, with field-perpendicular cross section 
L\ = (2-7r) 2 and field-parallel box size L\\. The choice 
of a rectangular box, as discussed in [l6| , allows for the 
excitation of elongated modes at large scales, necessary 
to avoid a long transition region from the forcing to the 
inertial interval, which can lead to inaccurate measure- 
ments of the spectral index. 

In order to achieve a steady state, the random forc- 
ing f is applied in Fourier space at wave-numbers 1 < 
k± < 2, (27r/X|i) < fey < (27r/Z/||)ni|, where nil controls 
the width of the force spectrum in ku . The Fourier coef- 
ficients inside that range are Gaussian random numbers 
with amplitude chosen so that the resulting rms velocity 
fluctuations are of order unity. The individual random 
values are refreshed independently for each mode on av- 
erage every r = 0.1 L±/v rms . As shown in [16|], the 
width of the field-parallel spectrum controls the critical 
balance at the forcing scale, and determines whether the 
turbulence is weak or strong. We define the Reynolds 
number as Re = (L±/2Tt)v rms /is. 

In the present simulations, we also introduce correla- 
tion between v and b, to investigate the role of cross- 
helicity. Such correlation is introduced at the forcing 
scales by controlling the correlation between the veloc- 
ity and magnetic field forces, f„ and f^. This is achieved 
by taking f as uncorrelated Gaussian random forces, 
-f ± + vV 2 z ± , (2) so that f„ = i (f+ + f"), f 6 = \ (f+ - f"), from which 



in which force and dissipation terms have been added to 
address the case of steadily driven turbulence, and we 
assume that viscosity is equal to resistivity. This set of 
equations is known as the Reduced MHD model (RMHD) 



it immediately follows that cross-helicity is controlled by 
the difference in the variances: (f„ ■ ft) = \ (a 2 + — crt), 
where o± = (If*! 2 ). It is convenient to define the pa- 
rameter a and the normalized cross helicity: 



26l |27I | . It is worth mentioning that the RMHD model 



was derived as an approximation of the full MHD equa- 
tions in the limit fey <C kx, and therefore it is applicable 
to strong turbulence. It has been recently realized that 
when used in a broader fcii — fej_ domain the system ([2]) 
describes the universal regime of weak Alfvenic turbu- 
lence [16( . This opens the possibility of effective analytic 
and numerical study of both weak and strong MHD tur- 
bulence in the same framework. In the present paper, we 



(o-l-a 2 _)/(a 2 + + a 2 _) = 2(i v -i b ), 


(3) 


H c /E=(E+-E-)/(E++E-) 


(4) 



Results. — Table U shows the summary of three repre- 
sentative simulations. Runs A, B, C were carried out 
with narrow fey -band forcing that produces critically bal- 
anced and strongly interacting large-scale modes. The 
energy spectra are shown in Fig. [TJ compensated by 
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FIG. 1: Spectra of balanced and imbalanced strong turbu- 
lence corresponding to runs A, B, CI (blue) and C2 (black) 
(see the text). In plots (B) and (C), the spectrum of z~ is 
arbitrarily offset in the vertical direction. 
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fcj_ . In run A with a = 0, we observe a balanced tur- 
bulence with the scaling close to E + ~ i? - ~ kj_ . 
When cross- helicity is introduced in Run B, we observe 
a slight steepening of E + ~ kT 1 ' 6 and a slight flatten- 
ing of E~ ~ fcj 1,35 . This behavior is justified as follows: 
since the total energy E = E + + E~ is kept constant, 
when the cross-helicity increases, the amplitude of z + in- 
creases at the expense of z~. Therefore, the nonlinear 
interaction of z + with z becomes weaker, resulting in a 
steepening of the spectrum. This steepening is however 
an artifact of a not optimal numerical setting. To sim- 
ulate this interaction correctly, we need to elongate the 
box in field-parallel direction so as to fit the eddies with 
longer parallel wavelengths at the forcing scales. As a re- 
suit, the E ± spectra get closer to k ± ' ; this is evident in 
Runs CI and C2. Note that the limit of very large cross- 
hclicity would require extremely long simulation box in 
order to observe the universal scaling behavior k ± . 

Discussion. — In this section we propose an explanation 
for the observed spectra. Our explanation essentially 
relies on the phenomenon of scale-dependent dynamic 
alignment. To understand how the alignment affects the 
energy spectrum, consider the eddy shown in Fig. [2] In 
this eddy fluctuations are aligned within the small an- 
gle 9\, while their directions and magnitudes change in 
an almost perpendicular direction. In the case of strong 
balanced turbulence, the nonlinear interaction in such an 
eddy is then reduced by a factor 9\ for both z + and z~ 




FIG. 2: A correlated region of (counter-) aligned magnetic 
and velocity fluctuations (red and blue vectors) at scale 
A = I/x/12, in a plane perpendicular to the strong guide 
field, in run A. The fluctuations are aligned predominantly in 
the x direction while their directions and amplitudes change 
predominantly in the y direction. 




FIG. 3: Sketch of dynamic alignment of magnetic and velocity 
fluctuations in a turbulent eddy. 



fields, and the corresponding nonlinear interaction time 
is estimated as t\ ~ 1/( z a ' kj_) ~ 1/( z a k±9\). The 
scaling of the fluctuating fields is then found from the re- 
quirement of constant energy fluxes: (z x ) 2 /t\ = const. 



One can argue 1J, |l5|, ll|| that the alignment angle de- 
creases with scale as 9\ oc A 1 ' 4 , in which case the field- 
perpendicular energy spectrum is E(k±) oc k J . 

In the case of strong imbalanced turbulence, the align- 
ment is still preserved. However, since the fields ampli- 
tudes are essentially different the alignment angles are 
different as well; we denote them 9t and 97, see Fig. [3] 
The assumption of the dynamic alignment then leads to 
the important geometric constraint: 9tzt ~ 97 z7, as 
is clear from Fig. [3J The depletion of nonlinear interac- 
tion is therefore different for z + and z~ fields, however, 
their nonlinear interaction times, r^ ~ V( z a ^±9 x ), are 
The requirement of constant energy fluxes 



const then ensures that z^ / z 



A/^A 



the same. 

y / e+/e , so both fields should have the same scaling, 
although different amplitudes. The geometric constraint 
then leads to 9^/97 ~ x /e~/e+, so the alignment angles 
should have the same scaling as well. 

To conclude this section we compare our results with 
recent analytic predictions of Lithwick et al. [20| , Beres- 
nyak and Lazarian 22j, and Chandran 2JJ. The main 
difference of our model with previous studies is that we 
include the phenomenon of dynamic alignment. In [20] 



the energy cascade times were assumed to be essen- 
tially different and the derived spectra had the form 



while our numerical results 
: 2 In [HI, EU, it was as- 



E+(k±) ex El oc k ± 5/3 , 

are more consistent with k 
sumed that the z + field undergoes a weak cascade, while 
z~ a strong cascade, leading to different spectra of z + 
and z - , which seems to be supported by numerical sim- 
ulations in [221] . We however note that the steepening 
of the z + and the flattening of the z~ spectra might be 
due to the high level of cross-helicity in these simula- 
tions (around p c ~ 0.97), which requires an extremely 
elongated simulation box in the field-parallel direction, 
cf. our Fig. QJB). In addition, the simulations in [22J are 
performed on lower resolution with hyperviscosity, which 
might alter the spectra due to a bottleneck effect. 




FIG. 4: Cosine of the alignment angle between v^ and b^ 
fluctuations in the guide-field perpendicular plane at scales 
A = L ± /6 (left), and A = L±/12 (right) in Run A. The right 
frame corresponds to the region delimited by the white square 
on the left side. 

Conclusion. — We have presented the results of numer- 
ical simulations of strong MHD turbulence in both bal- 
anced and imbalanced regimes. In the imbalanced tur- 
bulence, say with positive cross-helicity, the total energy 
spectrum E = E + +E~ is dominated by E + . Simulations 
in this case show a universal inertial-range regime: al- 
though the E + and E~ spectra have different amplitudes, 
their scaling is the same E + (k±) oc E~(k±) oc k ± . 

In the balanced turbulence, both spectra have the same 
amplitudes and scaling E + (k^) ~ E~(k±) ~ h]_ ' . 
This is consistent with the view that overall balanced 
turbulence can be imbalanced locally, creating patches 
(eddies) of positive and negative cross-helicity. In each 
of these regions the picture of imbalanced turbulence pre- 
sented above applies. When averaged over all the regions, 
the spectra of balanced turbulence are reproduced. 

The presented picture of MHD turbulence is consis- 
tent with the phenomenon of scale-dependent dynamic 
alignment 17J, 1 181 ] . which provides a natural explana- 
tion for the observed spectra. In this phenomenology, 
the configuration space splits into eddies with highly 
aligned and anti-aligned magnetic and velocity fluctua- 
tions, where nonlinear interactions are reduced. Fig. 0] 
shows the cosine of the alignment angle between veloc- 



ity and magnetic fluctuations for the balanced simula- 
tions (run A), at scales A = L±/6 and A = L±/12, The 
alignment possesses a hierarchical structure: inside small 
eddies there exist smaller and more anisotropic eddies 
(of both polarities), and so on. This hierarchical struc- 
ture, first observed by Mason & Cattaneo [unpublished, 
2006] , presents an interesting example of magnetic self- 
organization in a driven turbulent system. 
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